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ABSTRACT 


In this thesis, a performance analysis and improvement of a phase sampling inter- 
ferometer antenna system based on the Robust Symmetrical Number System (RSNS) in 
the presence of noise is investigated. Previous works have shown that the RSNS-based 
DF technique can provide high bearing resolution with a minimum number of antenna 
elements. However, the previous experimental data showed significant deviation from the 
theoretical results expected due to imperfections, errors, and noise. Therefore, an additive 
Gaussian noise model of RSNS-based DF was established and simulated. Simulation re- 
sults show that the presence of noise distorts the signal amplitudes used in the RSNS 
processor and causes degradation of the angle-of-arrival estimates. A performance analy- 
sis was undertaken by first introducing the quadrature modulation configuration into 
RSNS-based DF system, which provided a digital antenna approach for more flexibility 
in the signal processing. With a digital approach, variable resolution signal preprocessing 
can be employed, using a virtual channel concept. The virtual channel concept changes 
moduli values without changing the actual physical antenna element spacing. This attrac- 
tive property allows the RSNS algorithm to be implemented into existing antenna arrays 
and only requires modifying the antenna signal processor. Computer simulation results 
showed that the proposed method can successfully improve the system performance and 


also mitigate the effects of noise. 
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EXECUTIVE SUMMARY 


Radio direction finding (RDF) systems acquire the angle of arrival (AOA) of an 
incident electromagnetic wave (EM) by extracting the phase or time of arrival informa- 
tion from the intercepted waveform. The AOA data is prerequisite information in many 
modern technical applications, especially in crisis reaction areas. One of the more signifi- 
cant examples in the military sector is Early Warning Threat Detection, [1] which refers 
to methods that detect and identify the presence of an imminent threat of attack and im- 
mediately take countermeasures. The design of a RDF device may be divided into three 
categories: amplitude comparison, phase delay, and time delay. The phase delay tech- 
nique usually can generate an AOA resolution of +1 degree, which satisfies modern elec- 
tronic warfare (EW) requirements. Thus, most radio direction finding systems are based 


on phase-delay techniques. 


The direction finder based on the phase delay technique increases element spacing 
to improve the system’s resolution; however, it poses an ambiguity problem. It exists 
when a direction finder is composed of only two antenna elements whose separation is 
greater than half of the EM wavelength. There exists a variety of schemes to overcome 
the problem, such as The Maximum Likelihood (ML) Method, Eigenspace Methods, and 
Estimate of Signal Parameters via Rotational Invariance Technique (ESPRIT). However, 
due to the corresponding computational complexity, few practical systems have adopted 


these high resolution techniques. 


This thesis studied a new phase-sampled linear interferometer technique that can 
be easily implemented to provide a high-resolution, small baseline array. The approach 
relies on the unique properties of the robust symmetrical number system (RSNS) [7-10], 
which decomposes the spatial filtering operation into a number of parallel sub-operations 
(moduli) that are of small computational complexity. Each sub-operation is considered as 
a pre-configured interferometer that symmetrically folds the detected phase difference 


with folding period equal to 2Nm, where i=1,2,...,N, with N the number of interfer- 


XV 


ometers that are used within the linear array, and m, is the corresponding modulus of in- 


terferometer i. 


Previous experimental results from RSNS-based direction finders have shown the 
basic features of this signal processing technique. Since the previous designs were under 
noise-free theoretical derivations, the performances were marginal with some large AOA 
reporting errors. Therefore, the existing equations were reexamined and reconstructed 
from the standpoint of an additive Gaussian noise model. The system performance as- 
sessment for variable SNR was investigated and some fundamental relationships between 


SNR and the AOA resolution are presented. 


This study focused on developing a preprocessing step between the elements and 
the original RSNS processing circuit to improve the performance of the previous RSNS- 
based direction finder. Periodic element spacing in RSNS algorithm concept was investi- 
gated, which has the benefit to be implemented into existing antenna arrays with simply 


modifying the antenna signal processor. 


Finally, a variable resolution signal processing approach based on the virtual 
channel concept was introduced that provides an increased resolution when the SNR is 
high. Simulation results corroborated the analysis and clearly demonstrated the signifi- 
cant reduction in AOA estimation error that can be achieved with this new design con- 


cept. 


XVi 


I. INTRODUCTION 


A. RADIO DIRECTION FINDING SYSTEM 

Radio direction finding (RDF) systems acquire the angle of arrival (AOA) of an 
incident electromagnetic (EM) wave by extracting the phase or time of arrival informa- 
tion from the intercepted waveform. The AOA data is prerequisite information in many 
modern technical applications, especially in crisis reaction areas. One of the more signifi- 
cant examples in the military sector is Early Warning Threat Detection, [1] which refers 
to methods that detect and identify the presence of an imminent threat of attack and im- 
mediately take countermeasures. Therefore, numerous techniques have been developed to 


satisfy a wide variety of requirements. 


The design of a RDF device may be divided into three categories, depending on 


the fundamental measurement principles [2]: 
Category I: Amplitude comparison. 
Category II: Phase delay. 

Category III: Time delay. 


A Category I system converts the amplitude response of the antenna elements into 
voltages in order to obtain the AOA information. Figure 1 depicts the two types of an- 
tenna amplitude techniques used to obtain AOA information. They are referred to as the 
direct and comparative amplitude techniques. Here, @ represents the AOA to be meas- 


ured. 


For Categories II and III, at least two separate antenna elements with predeter- 
mined space between them are required. Figure 2 depicts the basic configuration for the 
phase-delay and time-delay techniques. We will assume a single frequency time- 
harmonic waveform. Since the electromagnetic wave experiences different propagation 
distances from an emitter to the individual antenna elements (except in the broadside 
case), it induces different phase terms at the antenna terminals. Equivalently, the wave 


arrives at each antenna element at a different time. The phase difference or time differ- 


ence is dependent on the element spacing, wavelength, and AOA. Since the first two fac- 


tors are known, the system can compute the AOA information. 


yY 





Horizontal E-field Dipole 
(a) (b) 


Figure 1. (a) Direct and (b) Comparative Amplitude Techniques. (After [2].) 
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Figure 2. Basic Phase-Delay/Time-Delay Techniques. (After [2].) 


Phase comparison systems usually can generate an AOA resolution of +1 degree, 
which satisfies modern electronic warfare (EW) requirements [1]. Thus, most radio direc- 
tion finding systems are based on phase-delay techniques. However, to improve the sys- 
tem’s resolution, the element spacing must be increased, which poses an inherent prob- 
lems of direction ambiguities. They exist for the phase-delay techniques when an inter- 
ferometer is composed of only two antenna elements whose separation is greater than 


half of the signal’s wavelength. These ambiguities cause error estimates in the angle of 
2 


arrival. There exists a variety of schemes to resolve the ambiguities, several of which are 


reviewed in the following paragraphs. 


The Maximum Likelihood (ML) Method [3] estimates the AOA from a given set 
of sampled phased data by maximizing the log-likelihood function. The likelihood func- 
tion is the joint probability density function of the sampled data, given the AOA, and 
viewed as a function of the desired variables (the AOA in this case). The ML methods 
yield nearly-optimal performance. Even when the signal-to-noise ratio (SNR) and the 
number of samples are small, they still have better performance compared to other meth- 


ods. ML methods, however, require significant computational resources. 


Eigenspace Methods rely on a property of the array covariance matrix. The space 
spanned by its eigenvectors may be partitioned in two subspaces, namely, the signal sub- 
space and the noise subspace when sources are incoherent. The steering vectors corre- 
sponding to the directional emitters are orthogonal to the noise subspace. In principle the 
Eigenspace methods search for directions such that the steering vectors associated with 
these directions are orthogonal to the noise subspace and are contained in the signal sub- 
space. The Multiple Signal Classification (MUSIC) [4] algorithm is an eigenspace 
method. Many improvements have been proposed to achieve better estimation perform- 
ance, such as Root-MUSIC [5]. The Root-MUSIC algorithm solves for the root of a 
polynomial instead of performing the conventional MUSIC search, and exhibits better 


performance and a more reasonable computational load. 


A third method, Estimate of Signal Parameters via Rotational Invariance Tech- 
nique (ESPRIT) [6] has a lower computation and storage requirement than other eigen- 
space techniques such as MUSIC. This advantage is achieved by using two identical ar- 
rays, in the sense that array elements need to form matched pairs with identical displace- 
ment by the same distance and in the same direction relative to the first element. An array 
geometry must be designed to have this property. Due to the corresponding computa- 


tional complexity, few practical systems have adopted these high resolution techniques. 


This research studied a new phase sampled linear interferometer technique that 
can be easily implemented to provide a high-resolution, small baseline array. The ap- 


proach relies on the unique properties of the robust symmetrical number systems (RSNS), 


3 


[7-10] which decompose the spatial filtering operation into a number of parallel sub- 
operations (moduli) that are of small computational complexity. Each sub-operation is 
considered as a pre-configured interferometer that symmetrically folds the detected phase 


difference with folding period equal to 2Nm, where i =1,2,...,N, with N the number of 
interferometers that are used within the linear array and m, is the corresponding moduli 


of interferometer i. 


The RSNS technique can be applied to either the amplitude or phase at the ele- 
ment output. The parallel use of phase waveforms increases the antenna resolution with- 
out increasing the folding rate of the system. Small amounts of analog-to-digital process- 
ing are required to quantize each folded phase response so that each sub-operation only 
requires a precision in accordance with its modulus. This feature permits the system to 
overcome phase errors induced by the array environment, including mutual coupling. A 
much higher AOA resolution is achieved by using more moduli in a system. Previous ex- 
perimental results from RSNS-based direction finder [7-10] have shown the effectiveness 
of this signal processing technique. The previous research presented a “mixed signal” 
processing approach. This thesis continued an examination of the performance of this 
novel signal processing under noise environments, and how a digital RSNS architecture 


might be implemented. 


B. PRINCIPAL CONTRIBUTIONS 

This thesis research focused on developing a preprocessing step between the ele- 
ments and the original RSNS processing circuit to improve the performance of the previ- 
ous RSNS-based direction finder. A new algorithm was established and new features 
were introduced, including the possibility of obtaining variable resolution with a fixed 


element spacing. 


The previous RSNS direction finders were designed under noise-free theoretical 
derivations and, therefore, the performances were marginal with some large AOA report- 
ing errors [7-10]. The existing equations were re-examined and re-constructed from the 
standpoint of an additive Gaussian (thermal) noise model. Computer simulations with 


noise were conducted in order to evaluate the previous system designs. The system per- 
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formance assessment for variable SNR was investigated. Some fundamental relationships 


between SNR and the AOA resolution are presented. 


During the course of the simulations, it was discovered that different moduli pairs 
under the same SNR condition have different performances. Since the selection of moduli 
affects the system performance, the possibility of changing moduli values without chang- 


ing the actual physical antenna element spacing was investigated. 


The possibility of employing equal spacing and still using a RSNS processing al- 
gorithm was also examined. Periodic element spacing is a very desirable characteristic. If 
the RSNS algorithm can be applied to periodic elements, then the concept can be imple- 
mented using any conventional array. It would simply be a matter of modifying the an- 


tenna signal processor. 


Finally, a variable resolution signal processing approach based on virtual channel 
concept is introduced that provides an increased resolution when the SNR is high. Simu- 
lation results corroborate the analysis and clearly demonstrated the significant reduction 


in AOA estimation error that can be achieved with this new design concept. 


Cc. THESIS OUTLINE 

The organization of this thesis is as follows. Chapter IT reviews the basic operat- 
ing principles of RDF systems, in general, and the two-element phase sampled interfer- 
ometer, in particular. The equations governing the sampling of an incident wave front are 


presented and the variables that affect the output are discussed. 


Chapter III is concerned with the RSNS theory. The application of RSNS theory 
to phase sampled interferometry is described, and the previous experimental results are 


presented. 


Chapter IV addresses the establishment of a noise model and the analysis of noise 
effects in a RSNS-based interferometer. The quadrature demodulation approach was in- 
vestigated to provide more flexibility on performance improvement. The possibility of 
applying the RSNS theory to uniformly spaced arrays was also studied. A variable resolu- 
tion signal processing method was proposed that provides a performance improvement in 


the presence of noise. 


Chapter V presents several case studies and computer simulations using the equa- 


tions derived in Chapter IV. Simulation results are discussed and analyzed. 


Chapter VI gives a summary of this study and recommendations for further re- 


search. 


The Appendix presents the MATLAB codes used to design, simulate, and evalu- 
ate the RSNS-based DF system. 


Il. PHASE DELAY-BASED INTERFEROMETRY 


In this chapter, we review the basic operating principles of RDF systems, in gen- 
eral, and the two-element phase sampled interferometer, in particular. The equations gov- 
erning the sampling of an incident wave front are presented and the variables that affect 


the output are discussed. 


ies THEORY 

At the outset, it is important to understand that under ideal circumstances, the 
electromagnetic wave incident on DF antennas exhibits a locally plane wave structure 
with linear polarization when the receiver is far enough from the transmitter and the an- 
tenna interactions with the platform are controlled. An ideal direction finder has the ca- 
pability of detecting signals over three-dimensional space. It is achieved by combining 
any two orthogonal, one-dimensional linear interferometers in a plane. Therefore, a brief 
description of the underlying theory of a one-dimensional linear interferometer is helpful 


in understanding the material in later chapters. 


A simple two-element linear interferometer based on phase delay interferometry 
is shown in Figure 3. Two identical antenna elements are separated by a distance d and 


the plane wave impinges at an incident angle 6 measured from the perpendicular to the 


baseline axis [7]. The angle @ takes on values from —90° to 90°. 


Incident Plane all 


a 
dsin(@) |-7 
a 


‘ 


Antenna 2 gt Re , Antenna 1 





Figure 3. Two-Element Interferometry. (After [7].) 
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Referring to Figure 3, the incident wave arrives first at antenna element | and then at 


element 2 after propagating an additional distance of dsin(@). The signals received by 


the two antenna elements have the following form 
V(@,t)=Vcos(a.tt+y,+¢,) (2.1) 


where i represents the antenna index number, V is signal amplitude response, q@, is the 
carrier frequency in radians/second, @¢, is the phase term relative to the origin, and y, is 


the phase term caused by the transmission line. The received signals are fed into a phase 
detector comprised of a mixer and a low-pass filter (LPF). After the phase detector mix- 


ing process, the received signal has the form 
Vv? 
Viv (Ot) =V (45t)V (Gt) =| cos (201+ Ly +X9)+cos(Ay+ Ag) | (2.2) 


where >) and A correspond to the sum and difference of the phase terms, respectively. 
After LPF filtering, the output signal from the phase detector has a DC voltage of 
the form 
Vv? 
V’ (d= 75 OS (Ay+A@) (2.3) 
where Ay comes from the transmission line length differences. The term A@ represents 


the phase differential between the two received signals, expressed as 


Ad=9,-, =kd sin (0) (2.4) 
where k is the wavenumber 
20 
k =— 2.5 
7 (2.5) 


and A is the wavelength at carrier frequency @.. 
It is reasonable to assume that Ay in Equation (2.3) can be measured and then 
eliminated. Let the term V* y 2 be normalized to unit amplitude, so that we obtain the fol- 


lowing final simple output signal from the mixing process: 


V..n (8) = cos (A@) = cos(kd sin()). (2.6) 


out 


B. AMBIGUITY AND ARTICULATION 
From Equation (2.4), the phase differential, Ag, repeats itself with a period 27. 


Thus the phase measurement becomes ambiguous at high frequencies or large spacings. 
Therefore, a single phase differential value corresponds to many angles of arrival. Ambi- 
guities are present when the distance between antenna elements is greater than half of the 


frequency wavelength, or d > 4/2. 
We take two representative numerical examples to demonstrate this criterion. 


First, consider the case of d = A/2, which gives Ag = sin(@). The phase differential 


versus angle is shown in Figure 4. There is a one-to-one mapping between AOA and 


phase difference. There are no AOA ambiguities created by spacing. 
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Figure 4. Phase Difference versus AOA for d = 4/2. (After [8].) 


Next, consider the case of d = A, which gives Ag = 27 sin(@). The phase differ- 


ence versus AOA, shown in Figure 5, suffers from ambiguities. The circles shown in 
Figure 5 have the same phase differential, which leads to the ambiguous estimation of 


angle of arrival. 
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Figure 5. Phase Difference versus AOA for d=. 


We can define the AOA measurement precision as the rate of change of the phase 


differential versus AOA, or 0 (AQ) / 00. The unit is degrees per degree [2]. Taking the 


derivative of both sides of Equation (2.4) yields 


0(Ag) | 
pe kd cos(@) (2.7) 


which shows that the measurement precision is proportional to the ratio of d// for con- 


stant carrier frequency and incident angle. Although a long baseline suffers from phase 


ambiguities, it increases the resolution of the AOA estimation. 


In practice, an interferometer DF system consists of dual baseline structures, so 
that the shorter baseline is used to overcome the ambiguities and the longer baseline is 


used to achieve high resolution [1]. 


C. FOLDING WAVEFORMS 
An examination of Equation (2.6) shows that the output voltage contains the plane 
wave AOA information in the form of a symmetrical folding waveform. This feature is 


illustrated in Figure 6 for the two representative spacing cases which have been dis- 


cussed. The number of folding periods n within the visible region (—2/2 <0 < 2/2) is 
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n=d/(A/2)=2d/A. (2.8) 
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Figure 6. Normalized Output versus Incident Angle for (a) d = 4/2 (b) d = A.(After 
[7].) 


By relating the element distances by a successive factor of two as shown in Figure 
7 and, appropriately shifting the output folding waveform of each baseline, it is possible 


to achieve a binary-resolution interferometer system [1]. For example, a 4-bit linear base- 
line interferometer consists of 5 antenna elements with spacings of 27A,2°A,2'A, and 
2° A. This kind of harmonic binary-related scheme is limited in its application due to re- 


quirements of a long baseline spacing to achieve high resolution. 


Antl Ant2 Ant3 Ant4 Ant$5 





Figure 7. Four-Bit Linear Baseline Interferometer. (After [1].) 


In summary, this chapter has discussed some basic direction finding techniques 
and their limitations. A signal preprocessing method employing the robust symmetrical 
number system is introduced in the next chapter. It can achieve a high resolution without 


large element spacing. 
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HI. ROBUST SYMMETRICAL NUMBER SYSTEM-BASED 
DIRECTION FINDING 


Since it is desirable to obtain high resolution AOA estimates without a long base- 
line, a signal preprocessing method employing the robust symmetrical number system 
(RSNS) was introduced in [7-10]. In this chapter, the RSNS theory is reviewed and pre- 


vious experimental data are presented and discussed. 


A. THE ROBUST SYMMETRICAL NUMBER SYSTEM 

The concept of RSNS-based preprocessing is that a robust symmetrical number 
sequence with proper shift value is used to represent a symmetrical folding waveform out 
of the antenna. Consider an N-channel RSNS, where N >1.The basic RSNS sequence has 


the following form: 


yeas | Pree 60s Brees Dereon (, aes (, ererreers Careeg (3.1) 
1 a a ae \ y 


where X,, represents the i” channel sequence (1 <i < N) with a corresponding modulus 


m,. Equation (3.1) originates from repeating each value N times in the following se- 


quence: 
[0,1, 2,....m—1,m,m—l,...,2,1] (3.2) 
where mm is a positive integer (m2 1). 
It has been shown that the period of a channel with modulus m, is [9] 


Pesys = 2Nm, 


L 


(3.3) 
and the fundamental period of an N-channel RSNS with modulus set {m,,m,, Bree my } is 
PFasys = 2N[m,,mM,,..., My | (3.4) 


where [m,,m,,...,My ] represents the least common multiple of m,,m,,...,m,. 


13 


The discrete states of each channel can expressed as [9] 








ae s,S1<Nm,+s,+1 

N 

aN) ON EN Sb 24 - 
oN | Nm, +5,+2<1<2Nm,+s,-1 


where g is the / state value of channel i, s, is the corresponding sequence shift defined 
as 


s, =0,1,2,---, N—I(mod NV) (3.6) 


and | x | indicates the greatest integer less than or equal to x. The sequence shift set 
Oe Sos a must be a complete residue system modulo N in order to have Gray-code 
properties. 

Table | lists discrete states within the RSNS for m, =3(s, =0) and m, =4 


(s, = 1). The column vectors, which consist of the values within each modulus, change 


one at a time between code positions (Gray-code property). 


ee ee 


srislelaalstel2[ali fore 
eb retstel ras helin| alisha) — 
Table 1. |Two-Channel RSNS for Modulus 3 and 4. 


It is desired to find the largest continuous set of column vectors shown in Table 1 








that are distinct. This sequence of distinct vectors forms the unambiguous output of the 
system and is defined as the system dynamic range, denoted by M. Studies [11] have 
shown that the selection of the sequence shift (s,) and permutations among the moduli 
has no effect on M, but may produce different start and stop points for the vectors corre- 
sponding to M. The dynamic range of the 2-channel example shown in Table 1 is 


M =15. That is, starting at vector (2,2) and ending at vector (3,1), all vectors are 
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unique. The symbol n re-indexes the sequence of RSNS vectors and counts the first vec- 


tor in the dynamic range as zero and the last vector in the dynamic range as M-1 or 14. 


The dynamic range of several particular moduli combinations have been formu- 


lated in closed form. For two channels, three combinations of two relatively prime 


moduli m, and m, can yield M as follows [10]: 
Case lI: m, 23 and m, =m, +1, then 


M =3(m,+m,)-6. (3.7) 


Case II: m, 25 and m, =m, +2, then 


M =3(m, +m,)—7. (3.8) 


Case III: m, 25 and m, =m,+C, and C 23, then 


M = 4m, + 2m, -2. (3.9) 
Checking the example shown in Table 1, which fits Equation (3.7), yields M =15. 


So far, we have reviewed the basic features of the RSNS. In the next section, the 


relationship between RSNS theory and DF systems is discussed. 


B. THE RSNS ANTENNA ARCHITECTURE 
There are three steps to integrate RSNS theory into a direction finding antenna 


design; these steps are presented in the following paragraphs. 


The first is to decide the physical antenna element spacing according to the se- 
lected moduli and operating frequency. Referring to Equation (3.3) and the definition of 


dynamic range, we define the number of folds within each modulus as 


n, =——=——. (3.10) 
Prsns 2Nm, 


Then combining Equation (2.8) and (3.10), we get the corresponding required antenna 


spacing for modulus m, 
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jeg es (3.11) 
2 4Nm, 


L 


In practice, the antenna element pattern degrades at wide angles so that the re- 
ceived waveforms are not ideal over a full 180 degrees. Therefore, the RSNS can be re- 
mapped to the reliable portion of the pattern to increase resolution [7]. Considering Equa- 


tion (2.6) with k and V_ held constant, we have 


dsin(@)=d’sin(6’) (3.12) 


where @ and d correspond to the maximum field of view (FOV) and the antenna spacing 
for the original (unscaled) configuration. The quantities 6’ and d’ are the maximum FOV 


and antenna spacing for the remapped (scaled) configuration. Rearranging Equation 


(3.12), the scale factor is defined as 





g-f = mel (3.13) 
d_ sin(@’) 
Combining Equations (3.11) and (3.13), we obtain the scaled antenna spacing 
; MA 
L=~2es 3.14 
(sdg=T—é G.14) 


Next, in an analog implementation, we need to design the comparators which en- 
code the normalized received waveforms shown in Equation (2.6) into RSNS form. The 


number of comparators required for each channel is equal to the corresponding modulus. 


It has been shown [7] that for a normalized folding waveform amplitude (-1<V,,, <1), 


the threshold for the j” comparator of modulus channel m, is 
SE 
Ves meos{ MI*P | 1< j<m,. (3.15) 


Whenever a folding waveform exceeds a comparator threshold, the comparator 
turns ON. A consequence of quantizing a folding waveform into discrete states is that the 
RSNS system displays any signal that falls within a discrete state bin as if it arrived at the 
bin center. Figure 8 illustrates discrete states of a two-channel RSNS with symmetrical 
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folding waveforms having moduli m, = 3(s, =0) and m, =4(s, =1). The thresholds 
Vx from Equation (3.15), shown on the vertical axis, encode the folding waveforms 
into the RSNS. The discrete states (integers) shown above the folding waveforms repre- 
sent the number of comparators in the ON state. 


|v-s —__—_—_—__} m=3,8;=0 
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Figure 8. RSNS Folding Waveforms for m, =3(s, =0) and m, =4(s, =1). (After 

[8].) 


Finally, we have to calculate the phase adjustment term for each channel in order 
to map the center of the dynamic range to the response of the antenna when the signal 
AOA is at broadside. Referring to Table 1, neither modulus has symmetry about the mid- 
point of the dynamic range, which is n=7. Therefore, a phase adjustment ¢ is added to 


Equation (2.6) giving 


V,,.(9) = cos(kd sin(@)+¢). (3.16) 


out 


When @=0, sin(@) =0, so V, 


out 


(9 =0) =cos(¢). The value of ¢ is selected so that 


V,, (@ =0) is encoded to the center bin of the system dynamic range. In the example 


out 


shown in Table 1, V. 


out 


(9 =0) should be encoded into 0 for m, and 3 for m,. 


At this point, we can obtain the relationship between the RSNS digital output and 
the angle of arrival. Using this relation, we further obtain the transfer function of the 
processor, which is a plot of true AOA versus estimated AOA. Studies have shown that 
the estimated AOA @’ for the u" bin is given as [7] 
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6’ =sin" Es -;| (3.17) 


and the AOA resolution r, for the u™ bin is given as 


Asin ee sin! eueaNl (3.18) 
M¢ M¢ 


In the next section, a design example is given to illustrate the important features 





of the RSNS processing method. 


C; DESIGN EXAMPLE AND EXPERIMENTAL RESULTS 

In this section, a two-channel, 6-bit RSNS antenna design originally given in [7] 
is described to illustrate the procedures in Section B. MATLAB programs are also devel- 
oped in order to provide essential computational modeling and simulation. All MATLAB 


programs are listed and explained in the Appendix. 


The schematic diagram shown in Figure 9 is for the design [13] of a 6-bit RSNS- 


based interferometer using m, =8(s,=0) and m, =17(s, =1). 
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Figure 9. Block Diagram of the RSNS-Based DF System. (From [13].) 
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The antenna elements are printed circuit dipoles operating at the frequency of 8.0 
GHz. For these RSNS parameters, we proceed to calculate the essential array configura- 


tion and transfer function with the assistance of MATLAB programs. 


1. Simulation Procedure 
Step 1. By Equation (3.9) we obtain the dynamic range, M =64. The RSNS vec- 


tor sequence is shown in Table 2. 
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Table 2. RSNS Vector Sequence Table for m, =8 and m, =17.(After [7].) 
Step 2. With a scale factor of & = 2/ V3, by Equation (3.13), the FOV changes 
from +90° to +60°. The corresponding element spacings by Equation (3.14) are 
d, =4.31 cmand d, = 2.03 cm. 


Step 3. By Equation (3.15) the corresponding thresholds of each channel’s com- 


parators are shown in Table 3. 
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Comparator | Channel m, | Channel m, | Comparator | Channel m, | Channel m, 


cos (337/34) 
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cos (1577/16) 


cos (1577/34) 





cos (1377/16) 


cos (1377/34) 





cos (117/16) 














\o 





Table 3. 








Comparator Threshold Table. (After [7].) 








Step 4. A phase adjuster ¢, for each channel is calculated using Equation (3.16). 


The normalized voltage output for channel m, after adding a phase adjustment is given as 


and for channel m, 


Vi (9) = cos (kd, sin(@) + 7/4) 
Viun2 (9) = cos (kd, sin(@)+ 7/2). 


(3.19) 


(3.20) 


Figure 10 depicts the folding output waveforms from Equations (3.19) and (3.20). 


The corresponding quantized folding output waveforms are plotted in Figure 11. 


Figure 10. 
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Folding Output Waveforms with Phase Adjustment. 
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Figure 11. Phase Adjusted Folding Waveforms after Comparators. 


Finally, by combining the information in Table 2 and Figure 11, and using Equa- 
tion (3.17), we can obtain a plot of estimated AOA versus true AOA as shown in Figure 
12. 
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Figure 12. Estimated AOA versus True AOA. 


Next, we compare the calculated results with measured experimental data and 


analyze the differences. 
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2. Experimental Results and Analysis 
A prototype circuit based on Figure 9 was built and tested [7, 8, 10]. The meas- 
ured AOA versus true AOA is shown in Figure 13. 


System Transfer Function 
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(degrees) 





Actual AOA (degrees) 
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Figure 13. Measured RSNS Transfer Function. (From [10].) 


Although we find that the transfer function shown in Figure 13 exhibits the basic 
features of the simulated (ideal) one shown in Figure 12, the performance is not good 
enough to achieve the DF objective. Possible error sources, such as comparators and mi- 
crowave components, have been discussed in [7, 8]. Here, we analyze the errors from a 
theoretical viewpoint. That is, an error analysis is performed on each parameter of the 
fundamental Equation (3.6) to show how errors effect the performance of the transfer 


function. 


First of all, we examine the antenna spacing d,. If the antenna spacing has a large 


fabrication tolerance (Ad =(d on )x100% ) or if the spacing does not follow 


theoretial 
Equation (3.11), it may cause errors similar to those shown in Figure 13. The error simu- 
lation result, shown in Figure 14, exhibits this phenomenon for the case of a spacing tol- 


erance Ad =+5%. 
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Figure 14. Transfer Function with Ad =+5%. 


The second parameter contributing to estimation errors is the phase adjuster c¢. 
Two cases of inaccurate phase adjustment were simulated. Figure 15 shows the errors 


resulting from an error tolerance of AC = (Cineoretiat — Sactuat )X 100% = 5%. Figure 16 shows 


the errors caused by mismatched ¢, which results in not mapping the broadside antenna 


response to the center of the dynamic range. 
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Figure 15. Transfer Function with Ac =+5%. 
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Figure 16. Transfer Function with Mismatched Phase Adjustment c. 


In this chapter, we have reviewed the basic features of the RSNS and the relation- 
ship between RSNS theory and DF system. An example of the previous RSNS-based DF 
system design demonstrated the design procedures. In the next chapter, we introduce the 
noise model into a RSNS-based DF system and examine the effects of noise on the trans- 


fer function. 
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IV. IMPROVEMENT IN RSNS DF SYSTEM DESIGN 


In this chapter, we develop several basic analytical theories that were used in the 
simulations of Chapter V. First, quadrature type demodulation is introduced, which pro- 
vides more flexibility in selection of a DF processing algorithm, and potential perform- 
ance improvement. Then, the additive Gaussian noise model is presented and integrated 
into the previous RSNS theory of Chapter II. For performance assessment, the Monte 
Carlo method was selected, which is widely used in modeling physical random processes. 
Finally, a virtual channel concept is introduced and employed to improve the RSNS- 
based DF system in the presence of noise. An overall RSNS system design block diagram 


is suggested at the conclusion of this chapter. 


A. QUADRATURE TYPE DEMODULATION 
1. Complex Representation of Band-pass Signal 


Bandpass signals are defined as signals whose spectrums are nonzero only within 


a region of + f, Hertz around a carrier frequency f,(=22/a,). Any bandpass signal can 


be expressed using complex notation as 


Xsana (Pat) = Re{x, (p)e’'} (4.1) 
where p represents the position information and x, (p) is a vector phasor that contains 


information on direction, magnitude, and phase. The phasor can be expressed as 





x, (p) =x, (p) + jXo(p)=|x, (p)|e*™ (4.2) 


where ¢() represents the phase term depending on the position. 
Using a trigonometric identity, we can rewrite Equation (4.1) as 
Xsana (Pst) = X, ( p) cos (@,t) — Xp ( p) sin (@,t) (4.3) 
where x, (p) is defined as the in-phase component and x, (p) is defined as the quadra- 
ture component. With respect to Equation (4.1), we call x, ( P) the baseband signal in 


complex form (or complex envelope), and e’’ the carrier signal in complex form. The 


25 


product of the two represents the modulated signal. The real part of this product is the 
transmitted signal. In many problems, it is convenient to suppress the time variation and 


represent the signal at any point by the vector phasor shown in Equation (4.2). 


2. Quadrature Type Demodulation 


As shown in Equation (4.2), the x, (p) and x, (p) terms preserve the phasor’s 


complete information. Therefore, quadrature type demodulation is widely used in many 
applications that need complete signal information. In order to further process antenna 
received waveforms, we employ the quadrature type demodulation structure shown in 


Figure 17 to replace the conventional structure (Figure 3). 


Assume that the signal received by antenna n is expressed as 


v(p,-t)=Re{A(p, eX} (4.4) 
where n=1,2,3 is the antenna element index and A( P,) is 
A(p,)=|A(p, Je” = 4, (P,) + Ao (P,)- (4.5) 


Figure 17 depicts a two-channel interferometer configuration using quadrature type de- 


modulation. 
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Figure 17. Quadrature Type Demodulation for a Two-Channel Interferometer. 
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A quadrature local oscillator signal (LO,) is obtained by shifting the local oscilla- 
tor signal, whose frequency is equal to the carrier frequency @., by 90 degrees as shown 
in Figure 17. That is 

Viog = C08(@,t+2/2)=—sin(@,t). (4.6) 
The received signal is fed to the mixer and mixed with the LO signal to get the in-phase 


term 


V1 (Pyst)=V(Pyst)cos(@,t) 





1 1 (4.7) 
=A; (p,)+5[ A (Px) €08(201) Ay (p, sin (20) 
and mixed with the quadrature signal to get the quadrature term 
Vo (p, .t)=—v(p,.t)sin(@,t) 
1 1 (4.8) 
= -5|4 (p,,)sin(2@,t)+ Ao (p, )eos(20,1) |+= Ag (p, )- 
After filtering with a low-pass filter (LPF), we obtain the baseband terms which are 
1 1 
(Px) = 54s (Px) = 5/A(Px)|e0s(0(P,))> (4.9) 
1 ik ; 
¥0 (Pn) =5 40 (P.)=5|A(P,)|sin(0(P, )). (4.10) 





2 


For a plane wave |A(p, ) 





=1 and therefore, we normalized Equations (4.9) and (4.10) to 


obtain the expressions 


Vn =V;(P,) =c0s(G(p, )), (4.11) 


Yon =Yo(P,) =sin(9(p, ))- oo) 
Equations (4.11) and (4.12) are referred to as the in-phase and quadrature terms of the 


signal received by antenna element n. 
Referring to Figure 17, the phase differential between antenna elements 1 and 3, 


A@,;, 1S expressed as 
AG; = 9(P,)— (D3) = kd, sin (8) (4.13) 


pai 


where @ is the incident angle measured from the perpendicular to the baseline (broad- 
side), and k is the wavenumber. The phase of antenna element 3 is established as the ref- 


erence phase, and d, is the distance between element 1 and the reference element 3. A 
similar procedure for antenna elements 2 and 3 produces 
Ad,, = 9( Pp. )—O(p,;) = kd, sin(@). (4.14) 


At this point, an in-phase output for a channel composed of antenna element n 


and reference antenna element r is defined as 
V,, (8) = cos (kd, sin(8)) (4.15) 
where d, is the distance between element n and the reference element r. Correspond- 


ingly, a quadrature output is defined as 


Vv 


on (9) = sin (kd, sin (8)). (4.16) 


Notice that Equations (4.15) and (4.16) can be related to Equations (4.11) and (4.12) as 


follows: 
V,, (8) = cos (kd, sin(8)) =cos(AQ,, ) 
=cos(¢(p,,))cos(@(p, ))+sin(¢(p, ))sin(¢(p, )) (4.17) 
= ViVi + VonVor 
and 


V,, (@) =sin (kd, sin(@)) = sin(AQ,, ) 
=sin(¢(p, ))cos(4(p, ))—cos((p, )) sin(¢(p,)) (4.18) 


- VonVir ~ VinVor : 


In a practical application, both terms for a channel can be computed in a channel proces- 
sor as shown in Figure 18. After obtaining I and Q for each channel, the processor output 


goes through the same processing as the conventional configuration does (Figure 3). 
Referring to Equation (3.16), a phase adjustment ¢, is added to channel 7 in order 


to map the channel output to the center bin of the system dynamic range when 0 = 0°. 


That is, 
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V, (8) =cos(kd, sin(@)+¢,). (4.19) 
Using a trigonometric identity, we can rewrite Equation (4.19) as follows: 
V, (@) = cos (kd, sin(@) +¢;) 
=cos(kd, sin(@))cos ¢, — sin (kd, sin (8) )sin ¢, (4.20) 
=V,,(@)cos¢, —V,,(@)sing,. 
Equation (4.20) shows that phase adjustments can be implemented in digital form instead 


of using a hardware phase shifter. That is, we can acquire the phase adjustment for each 


channel modulus in advance and save it in memory. Figure 18 shows the implementation 


based on Equation (4.20). 
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Figure 18. Quadrature Type Demodulation for the Second Channel Including a Digi- 
tal Phase Adjustment. 
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B. NOISE MODEL 

1. Additive Gaussian Noise Model 

The additive Gaussian noise assumption has been widely adopted for modeling 
the underlying random noise in many different applications, since it usually reduces the 
complexity of the problem from both theoretical and empirical standpoints. In this section, 
we add a bandpass random process in the processing algorithm to model the noise envi- 


ronment. 


Consider a bandpass noise signal in the following form: 


N(t)=Re{|A(r)|ee't (4.21) 





where N (t) is Gaussian noise, A(t) is a random process with Rayleigh first-order prob- 
ability density function (PDF), and g(t) is a random process with a uniform first-order 
PDF over the range (0, 27]. It is convenient to use the equivalent form: 


N(t)=N, (t)cos(@,t)— No (t)sin(@,t) (4.22) 


where N,(t) and N,(t) are given by 
N, (t) =|A(t)|cos(g(t)) (4.23) 


No (t) =|A(t)]sin(e(t)). (4.24) 
Combining Equations (4.7), (4.8), (4.23), and (4.24), we obtain a Gaussian noise model 
for the RSNS-based DF system shown in Figure 17. That is 


Vinoise (Past) =V( Pat +N; (t) (4.25) 


Vonoise (Puot) =V( Py ot) + No (t) (4.26) 


2: Monte Carlo Simulation Method 


Monte Carlo (MC) methods are based on the use of probability and statistics to 
model a random process by repeated trials. It was named by S. Ulam, after the city of 
Monaco, Monte Carlo, where the primary attractions are casinos providing games of 


chance. MC methods are applied in many cases where a closed-form statistical analysis is 
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mathematically complex or too difficult to perform. In order to adopt the MC method to 
estimate the performance of a RSNS-based DF system in the presence of noise and inter- 
ference, we first review the MC algorithm. The pertinent equations are derived and then 
implemented in MATLAB programs to obtain the estimation. All MATLAB programs 
are listed and explained in the Appendix. 

A random variable X that is used to describe a signal S$ distorted by noise W is 
expressed as 

xX =S+W (4.27) 

where W is a zero mean Gaussian random variable with variance 0” =1. Suppose that 


the probability of X >m fora given S' is expressed as 
P(m)=P(X >mlS) (4.28) 
where m is a constant. In order to estimate P(m), we generate a sequence of statistically 


independent and identically distributed Gaussian random variables w(n) with zero-mean 
and unit-variance, where n =1,2,3,...,N. Then, by adding a constant value s to w(n), we 


obtain a new set of random variables 
x(n)=st+w(n). (4.29) 


Assume that a new random variable y(n) is defined as 


7 0, ifx(n)<m 
vio)=fp if x(n)>m ; ee) 


Then, the estimation of P(m) is 


P(m)=—Y" y(n), (4.31) 


n=l 


For the precision of Equation (4.31), we need to calculate the mean and the vari- 


ance of P(m). The mean of P(m) is [14] 
A 1Q 1 
E| Pom) = LED ()] = 57 NP(m) = P(m) (4.32) 
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where E [e] denotes the expectation operation. The variance is 








A 2 
Fo) = Le (m) oe (m) | 
ee 1 (4.33) 
= E| P*(m)|- P? (m)= ay Pm)[1-P(m) 
Next, suppose that the precision of the estimation 1s expressed as 
P(|P(m)—P(m)|<e)2@ (4.34) 
where € is the relative error, defined as 
P(m)—P(m 
Psa) (4.35) 
P(m) 
and @ is the desired outcome percentage. By Chebyshev’s inequality, 
o 
P(|P(m)- P(m)]<e)21-—. (4.36) 
E 
Referring to Equation (4.34), it gives 
gia (4.37) 
Substituting Equation (4.33) into Equation (4.37), we obtain 
P 1-P 


(l-a)e* 
which is the required number of trials for Equation (4.34) to hold. 


Now we sum up the MC algorithm as follows: 


1. Set up initial number of trials, N, based on the desirable relative error, €, and per- 


centage, @. 
N 
2. Generate a sequence of random numbers { y (n)} 


n=l" 


3. Compute the number of { y (n)} that exceed m, denoted by N,. 


4. Estimate the probability P(m) = a, 


5. Use the estimate P (m) in Equation (4.38) to check if Equation (4.38) still holds. 


If yes, go to step 6; otherwise, calculate the corresponding N by Equation (4.38) 
then go back to step 2. 


6. Report P(m), then end the program. 


In the next chapter, this algorithm is used to assess the RSNS-based DF system’s 


probability of correctly estimating the signal’s incident angle in the presence of noise. 


Cc; VARIABLE RESOLUTION SIGNAL PROCESSING 
1, Equal Spacing Antenna Arrangement 


An array antenna is comprised of a collection of elements positioned in space and 
interconnected. By the way of electronic scanning, array antennas can achieve better per- 
formance than mechanical scanning antennas. For this reason, array antennas play an im- 


portant role in modern EM applications. 


The fundamental configuration of array antennas is the linear equispaced array in 
which radiation elements lie equispaced along a straight line as shown in Figure 19. Since 
this configuration is simple to implement with hardware, the linear equispaced array an- 
tennas are more widely used than any other array configuration. Therefore, it would be a 
major advantage if the RSNS processing technique could be applied to periodically 
spaced elements. If this were possible, then the RSNS method could be incorporated into 
existing antenna arrays, and there would be no requirement to design a custom array with 


unequally spaced elements. 





Figure 19. Linear Equispaced Array Antenna. 
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From Equation (3.11), when the frequency is constant, we obtain 


ia, (4.39) 


i m, 
Equally spaced elements require 


d, =id, (4.40) 


where i = 2,3,4,...,N for N channels. Combining Equations (4.39) and (4.40), we obtain 


xj, (4.41) 


Since the RSNS is based on a number of pairwise relatively prime (PRP) moduli 
[11], Equation (4.41) does not hold. The only way to satisfy Equation (4.41) is to extend 


the RSNS definition to non-PRP moduli. The cost is a decrease in the dynamic range of 


the RSNS. Table 4 shows a comparison of dynamic range values for both non-PRP and 


PRP RSNS. 


Non-PRP RSNS PRP RSNS 





























Table 4. |. Non-PRP versus PRP RSNS 


2. Virtual Channel 


Consider the two-channel interferometer with quadrature type demodulation 


shown in Figure 17. The moduli pair of m, and m, are selected for a two-channel RSNS 


whose dynamic range 1s M. As shown in Equation (4.17) and (4.18), the corresponding 


in-phase and quadrature terms for each channel are 


V,, =cos(kd, sin(8)), (4.42) 
Vo, = sin (kd, sin(8)), (4.43) 
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V,. =cos(kd, sin(@)), (4.44) 
Vo. = sin(kd, sin(@)). (4.45) 


By properly combining Equations (4.42) to (4.45), we obtain the following equations: 


Vii =ViVi2 —VoVo2 = c0s(k (d, +d, )sin(@)) =cos(kd/sin(@)), (4.46) 
Vis =V Vin +VoVo2 = c08(k(d, —d,)sin(@)) =cos(kd; sin(A)), (4.47) 
Vor =VoV 2 t+ViiVo2 = sin (k(d, +d, )sin(@)) = sin (kd/sin (8)), (4.48) 
Vo2 =VoV 2 —ViiVo2 = sin(k(d, —d,) sin(@)) = sin(kd; sin(0)). (4.49) 


Note that V; (i =1,2) can be regarded as the in-phase term of a virtual channel with ele- 
ments spaced d;. Similarly, Vj, is the quadrature term of virtual channel with elements 


spaced d/. 
Now suppose that 
mm, = he, (4.50) 


Then using Equations (3.14), (4.46), and (4.47): 





peed deen _AME 





dj =d,+d, = = 4.51 
Sie PON NCH “iy 22N m, $ 2 2Nm, en 
, AM({1 1 AMh-1, AME 

d, =d,—d, = c= g= cf (4.52) 

22N\m m, 22N m, 2 2Nm, 


which give the relationships for the virtual spacing to a new modular pair {m; ‘ m, , hav- 


ing the dynamic range M’, and the scale factor &’. Rearranging Equations (4.51) and 
(4.52) gives 
M(h+1), M’,, 
==; (4.53) 


mM, m, 
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M (h-1) 





=e. (4.54) 


2 


Then, dividing Equation (4.53) by Equation (4.54) 





tN, (4.55) 
h-1 m, 


These relationships allow RSNS processing to be performed for elements that do 
not meet the physical spacing requirements set forth in Chapter III. The overall modified 
processing design for this so-called variable resolution approach is described in the next 


section. 


D. OVERALL DESIGN 


From a system level view, a RDF device consists of four essential elements [2]. A 


block diagram of a generalized RDF device is shown in Figure 20. 






EM waves at AOA AOA Read-out Display 


a ee Receiver For Post-receiving Processor 
: RF/IF/Baseband Processing For AOA Information 


Figure 20. Functional Blocks for a Direction Finder. (After [2].) 





The first box is the antenna, which extracts electromagnetic energy and converts it 
into a signal containing the AOA information. The second box is the receiver, which 
converts, amplifies and processes the radio-frequency signal to intermediate frequency 
(IF) or baseband. The third box is a post-receiver processor, which further processes the 
signal in order to extract the AOA information. The fourth box is the display device that 


indicates the AOA. 


Note that the operations in Figure 17 and 18 are essentially the first two boxes of 
Figure 20. The RSNS decoding block in Figure 9 is equivalent to the third box of Figure 
20. It is at this point that we implement the virtual channel concept and integrate it into 
RSNS decoding block in order to improve performance without changing the array struc- 


ture. The flowchart of the variable resolution technique is shown in Figure 21. 
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: Receiver For’ : 
: RF/IF/Baseband : 
Processing: 






[AOA Read-out 
: Display 





: Post-receiving Processor For AOA information 


Figure 21. The Flowchart of the Variable Resolution Signal Processing Technique 


Referring to Equation (3.15), the higher m,, the more quantization levels are re- 


quired. Thus, the spacing between levels is smaller and therefore the coding operation is 
more susceptible to noise effects. However, I-Q detection allows for the possibility of a 
variable resolution processor. Initially we select low-resolution moduli so that the direc- 


tion finder is reliable in a low-SNR environment. This provides a rough estimate of the 


AOA 6. That is, the bin number of the AOA is known, but the bin width is large. Next 


we use high-resolution moduli to obtain a more accurate angle estimate 6, within the 


large bin width of the first estimate. We decide which estimated AOA is displayed, based 
on the reliability of its measurement. The second higher resolution estimate is only re- 
ported if its probability of being correct exceeds a predefined threshold. If the fine esti- 
mate angle is sufficiently reliable, then it is reported; otherwise the coarse one is reported. 
With this approach, under a low-SNR condition, we still obtain a coarse AOA estimate, 


but when the SNR is high a much more accurate AOA estimate is possible. 


In this chapter, several theoretical aspects of direction finding, signals and noise, 
and the basic RSNS DF processing have been discussed. A new variable resolution RSNS 
processor has been described. In the next chapter, simulations based on the techniques 
discussed in this chapter are used to predict the performance of RSNS direction finder 


under various operating parameters. 
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V. SIMULATION, RESULTS, AND ANALYSIS 


In this chapter, we present computer simulations to predict the performance of 
RSNS-based DF systems in the presence of noise. The approach based on the variable 
resolution signal processing technique is demonstrated and compared with the original 


RSNS design. 


A. NOISE MODEL FOR A RSNS-BASED DF SYSTEM 


In this section, we use the additive Gaussian noise model presented in Chapter [V 
to theoretically assess the performance of the interferometer shown in Figure 9. 
MATLAB programs using the Monte Carlo algorithm were developed to accomplish the 
assessment. All MATLAB programs are listed and briefly explained in Appendix. 


At the outset, we introduce the simulation parameters that are used throughout the 
whole chapter. The relative error € was set at 3%, and the desired percentage @ was set 
at 90%. Therefore, the required number of trials was calculated by Equation (4.38) and 
dependent on the corresponding value of probability. Noise was added to the folding 
waveforms out of each channel processor, which were in the range -1V <V,, <+1V and 
-1V< Vo = +1V, where i=1,2 is the channel index. The variance of the noise was de- 
termined from the specified SNR by 


Vv’ 
20 


SNR = 





(5.1) 


2 


where V, is channel output defined in Equation (4.20). 


After noise was added, the channel output voltage was quantized and the AOA 
bin computed. The estimated AOA bin was obtained by selecting the bin with the highest 
probability. The corresponding estimated AOA was obtained by converting the estimated 
bin to angle units using Equation (3.17). If the quantized output value could not be en- 
coded into any bins, no estimated AOA bin was obtained, and the estimated angle was set 
at zero. The probability of a correct bin was defined as the probability of the estimated 
bin being equal to the theoretical bin. The probability of a bin error was defined as one 


hundred percent minus the probability of the correct bin. 
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First, we examine the estimated AOA versus true AOA for moduli {8,17} with 
the SNR being set at 20 dB. Figure 22 shows the simulation results. Observing Figure 22, 
we find that the transfer function “fades” in two symmetrical regions about 0°. That is, 


under the same operating condition, some spatial regions may be estimated incorrectly 
with higher probability. This observation may be helpful in developing a solution to the 


“fading” problem. 
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Figure 22. Transfer Function for Moduli{8,17}, ¢ = a V3, and SNR = 20dB. 


Next, we investigated the probability of correct bin estimation over a FOV with 
variable SNR. This investigation allowed us to determine the RSNS processor operating 
environment. Figure 23 shows the simulation results. From this figure, we find that it re- 


quires approximately a SNR of 30 dB to achieve reliable estimation, which we define as 


the probability of choosing the correct bin being greater than 50%. 
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Figure 23. Probability of Correct Bin Estimation for Moduli {8,17} and & = 2//3 
with Variable SNR. 


Since the resolution of a RSNS DF only depends on two parameters: the dynamic 
range M and scale factor & [7-11], we simulated the DF processing for different moduli 


pairs which have the same M and & to see how noise affects their performance. Figures 


24 to 26 show the simulation results for three moduli pairs with SNR = 30 dB and 


E=2/V3. 
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Figure 24. Noise Model for Moduli 15; 23} (a) Probability of Bin Error and (b) Quan- 


tized Output for Modulus m, =5. 
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Figure 25. Noise Model for Moduli {8, 17} (a) Probability of Bin Error and (b) Quan- 
tized Output for Modulus m, =8. 





>+- MU oO fF A DW 
oO oO oO oO fo 
TTT TT 
jg 








Probability of Bin Error (%) 


-60 -40 -20 0 20 40 60 
(a) Angle of Arrival (degrees) 





fare 











Quantized Output 
O-]+-$NWHRUONDWOO 
T T T T T T T T T 








-60 -40 -20 0 20 40 60 
(b) Angle of Arrival (degrees) 
Figure 26. Noise Model for Moduli {10,13} (a) Probability of Bin Error and (b) Quan- 


tized Output for Modulus m, = 10. 


From the figures, we find that the probability of a bin error in all three cases is 
roughly the same, but the fading regions (1.e., areas of high probability of bin error) are 
quite different. The phenomenon is that the greater the difference is between the two 


moduli values, the closer together the fading regions appear. Furthermore, we observe 
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that the folding waveform for the smaller modulus in each case is correlated with the 
shape of the error curve. When the smaller modulus increases, this effect is reduced. 
Therefore, we conclude that different moduli with the same M , 6, and SNR have the 
same correct estimation probability but in different FOV regions. 

Not only is the probability of a correct bin important, but also the magnitude of 
the error. Figure 27 compares the simulation results for a noise-free RSNS processor and 
one with additive Gaussian noise (SNR = 20dB). Comparing Figure 27 with the 20 dB 
curve in Figure 23, it shows that although some regions have low probability of correct 
bin estimation, the processor still can get useful AOA estimations. Hence, this property 
may be helpful for improving the design. Referring to Figure 13, Figure 27 may explain 


why the experimental results cause huge errors at some angles. 


40 





30/7 


207 


Estimated Error (degrees) 
[= ) 














— Noise Added 
‘== Noise Free 
40 L 1 L 1 
-60 -40 -20 0 20 40 60 
Angle of Arrival (degrees) 





Figure 27. Estimated Error Comparison for Moduli {8,17}, ¢ = 2/ Neo 


B. VARIABLE RESOLUTION SIGNAL PROCESSING 
1. PRP RSNS 
Let us start with the PRP RSNS. We continue using the moduli m, =8, and 


m, =17 to illustrate the virtual channel concept and show how the concept increases the 


system resolution. We shall use a frequency, f, =8 GHz, scale factor, € =1, for which 


A 


the FOV is 90%, and dynamic range, M = 64. 
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By Equation (4.50), we obtain a modular ratio, h =17/8. Substituting the ratio h 
into Equation (4.55), we get the new moduli m; =9 and m’, = 25. The new dynamic range 
is acquired from Equation (3.9), which is M’ =84. Notice that it is higher than the dy- 
namic range of original moduli M. By Equation (4.53), the new scale factor is 
&’=150/119 which implies that the FOV is about +53’. Figure 28 shows the transfer 


function simulated for the new design settings. 
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Figure 28. Transfer Function for Moduli {9,25} and ¢’ =150/119. 


Although we could redesign the array for a new & equal to &’, we also must 
change the antenna element spacings. In the virtual channel approach, we not only in- 
crease resolution with a higher dynamic range value, but also keep the same physical 
element spacings. Therefore, if the signal arrives within the new FOV, the variable reso- 


lution approach can be used to reduce the estimation errors. 


At this point, we have shown the new design approach in the noise-free operating 
condition. We continued to investigate it under the noise-added operating condition. First, 
we examined the probability of correct bin estimation. Figure 29 shows a comparison of 
the two systems for a SNR of 30 dB. Notice that the two systems have slightly different 
fading regions (1.e., low probability regions) that can be used to compensate for each 


other’s fading region. 
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Figure 29. Comparison for Moduli{8,17} and Moduli {9, 25}. 


Next, we integrated the virtual channel concept into our incident angle estimation 


algorithm, whose flowchart is shown in Figure 21. Figure 30 shows a comparison be- 


tween the two methods. Clearly the virtual channel has resulted in smaller error than the 


original RSNS moduli. The SNR for this case is 30 dB, which is considered a relatively 
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We have shown that the virtual channel does improve our direction finding accu- 


racy for high SNR. However, as mentioned before, a low SNR may drop the probability 
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of correct bin estimation significantly (shown in Figure 23). That is, if the operating SNR 
is not high enough, then the virtual channel processing may not contribute much. Figure 
31 demonstrates this case. For a 20 dB SNR, the virtual RSNS processor has roughly the 


same performance as the original RSNS processor. 
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Figure 31. Comparison Between Original Method and Virtual Channel Method 


Thus, the original moduli determine the noise tolerance of the system. It is desir- 
able to have moduli providing a robust estimation in low SNR environments while hav- 


ing a high resolution capability. 


2. Non-PRP RSNS 
The example in Figure 31 belongs to Case III defined in Chapter II. We now ex- 


amine the other two cases to see whether the virtual channel processing can be applied. 


First we check Case I, which is m, 23 and m, =m, +1. Referring to Equation 


(4.55), we obtain 





Dee (5.2) 
m, 1 
For the Case II, which is m, 25 and m, =m, +2, we obtain 
m, _ 2m +2_ mth (5.3) 





m 2 1 
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Since both cases result in a modulus equal to one, we could not process the signal 
further, unlike Case III. In other words, the virtual channel concept cannot be applied. 
However, if we extend the RSNS definition to non-PRP moduli, then the virtual channel 


processing can be applied. Two numerical examples illuminate this fact. 
For Case I, moduli m, =5 and m, =6 with scale factor € =1 is considered. The 


corresponding dynamic range is M =27, calculated by Equation (3.7). By Equation (5.1), 


we obtain 





2 == (5.4) 


where the corresponding dynamic range is M’ =47. Referring to Equation (4.50), we ob- 


tain the new scale factor 


_ 99 


’— "7" e] 5.5 
oaae (5.5) 
which means FOV > +90. That is, not all RSNS bins lie in the visible region. For this 


case, the effective dynamic range is M’ = 21(=33-—13+1). Figure 32 shows the results. 
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Figure 32. RSNS Code versus AOA for Moduli {2,22}. 
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For Case II, moduli m, =5 and m, =7 with scale factor € =1 is considered. The 


corresponding dynamic range is M = 28, calculated by Equation (3.7). By Equation (5.1) 


we obtain 





(5.6) 


where the corresponding dynamic range is M’ =27. Referring to Equation (4.50), we ob- 
; yp a2 : ; 
tain the new scale factor €’ = B <1, which means the FOV > +90°. Figure 33 shows the 


results. 
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Figure 33. RSNS Code versus AOA for Moduli {2,12}. 


For both cases, although the virtual channel concept still holds, it does not con- 
tribute to a higher resolution AOA estimate. That is, the new dynamic ranges are lower 
than original RSNS dynamic range. In next subsection, we investigate the possibility of 
equal spacings and analyze the potential benefits. 


3. Equi-spaced RSNS 


As mentioned in Chapter IV, accepting non-PRP RSNS moduli allows equal ele- 


ment spacings, but at the expense of decreasing dynamic range. Thus, we try to integrate 
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the virtual channel processing into an equi-spaced RSNS array and observe whether we 


can obtain other benefits to compensate for the dynamic range decrease. 
Consider the following scenario: 
e Two channels (three elements), where the elements have equal spacing. The 
frequency is f, =8 GHz, and the SNR = 20 dB. 
e ARSNS design based on moduli m, =8 and m, =16 with a scale fac- 
tor € =1. The corresponding dynamic range is M =35. 


By Equation (4.47), we obtain h = 2. Substituting the ratio into Equation (4.56), 


we get the ratio of the new moduli as 





(5.7) 








, 
m —_ 1 





Table 5. Moduli for 


, 


Ms 


According to the variable resolution processing algorithm shown in Figure 21, ini- 
tially we select a low resolution moduli to get a coarse AOA estimate which is more ro- 
bust. Thus, the moduli m/ = 3 and m’, =9 are selected. The noise model performance 
comparison between the original moduli and the new moduli is shown in Figure 34. The 
data shows that the original moduli have reliable estimations only over a small region 
while the new moduli yield 50% or higher over most of the FOV. This result is expected; 
that is, the new moduli provide a better estimate than original one under a low SNR situa- 


tion. 
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Figure 34. Noise Model for Moduli {3,9} and {8,16}. 


Next, let us examine the scale factor, which affects the FOV. By Equation (4.53), 
we obtain €’ =105/112 which means FOV > +90°. Figure 35 shows the simulation result. 
In this case all RSNS bins are used, but parts of the first and last bins lie outside of the 


visible region. 
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Figure 35. RSNS Code versus AOA for Moduli {3,9}. 


Next we examine another case of interest. That is virtual moduli that yield a reso- 


lution higher than the original moduli. The moduli mj =8 and m’, = 24 are chosen, with a 
50 


SNR of 20 dB. The noise model performance comparison between the original moduli 


and the new moduli is shown in Figure 36. 
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Figure 36. Noise Model for Moduli {8,16} and {8,24}. 


Figure 36 again demonstrates the phenomenon that the RSNS processor with a 
higher dynamic range suffers from noise more easily than one with lower dynamic range. 


Notice that for moduli m; =8 and m’, = 24, the corresponding scale factor is €’ = 105/102 
which means FOV ~ +74. 

The above three RSNS designs share the same element spacing, but each one has 
a different level of resolution and robustness. We can employ all three moduli at the same 


time or only select one of them depending on the operating SNR without changing the 


physical element spacings. The moduli m/ = 3 and m’, =9 provide the most robust esti- 
mations at lowest SNR. The moduli m, =8 and m, =16 provide an intermediate level of 
estimation, while the moduli m/ =8 and m’, = 24 provide the most accurate estimation. 
With single element spacing, a three-level progressive angle strategy can be em- 
ployed following the flowchart of Figure 21. First, we determine the error versus FOV for 
the RSNS moduli m, =8 and m, =16, which is shown in Figure 37, and use it as the per- 


formance baseline. 
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Figure 37. Estimated Error for Moduli {8,16}. 


Then, following the algorithm of Figure 21, Level one estimation is conducted 
with moduli m; =3 and m’, =9. The simulation result is shown in Figure 38. Compared 
with the original moduli in Figure 37, Figure 38 shows the significant improvement be- 


tween 40° and 60° from broadside. 
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Figure 38. Comparison Between Original Method and Level One 


After Level one processing, the signal continues to be processed by a second set 


of moduli, m, =8 and m, =16, which is called Level two estimation. Figure 39 shows the 


a2 


result of this step. This time the improvement is not as dramatic, but there is still some 
improvement. 
The last level of processing (Level three) is done with moduli m; =8 and 
m,, = 24. Figure 40 shows the final level processing results. At this point, only the central 
regions are improved from Level two. Therefore, for this selected scenario, two level 
processing is enough. 
When the operating SNR is high enough, we can directly employ the moduli 
m, =8 and m, =16 at Level one processing. The drawbacks of adopting non-PRP RSNS 


coding are compensated for by this attractive variable resolution signal processing capa- 
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Figure 39. Estimated Error versus AOA for (a) Level One (b) Level Two 
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In this chapter, we have discussed the RSNS-based DF system in the presence of 
the additive Gaussian noise. Simulation results corroborated the Chapter IV analysis and 
clearly demonstrated the decrease in performance. Also, we have shown the effective re- 
duction in angle of arrival estimation error that can be achieved with variable resolution 
signal processing. In the next chapter, we summarize the results of this study and intro- 


duce some possible topics for future research. 
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CONCLUSIONS AND FUTURE WORK 


This research has analyzed the RSNS-based DF system from the viewpoint of an 
additive Gaussian noise model. Computer simulations showed that the presence of noise 
distorts the signal amplitudes applied to the RSNS encoding and caused a decrease in the 


AOA accuracy. Also the performances of different moduli with the same parameters have 
been assessed. Simulations concluded that different moduli with the same M ,¢, and SNR 


have the same correct estimation probabilities, but in different FOV regions. 


The possibility of employing equal element spacing in the RSNS processing algo- 
rithm was investigated by extending the RSNS definition from PRP to non-PRP. The 
drawback of adopting non-PRP RSNS coding is a decrease in dynamic range. However, 
it was compensated by the development of a variable resolution signal processing ap- 
proach. This approach leads to the virtual channel concept, which refers to changing the 
moduli values without changing the actual physical antenna element spacing. This attrac- 
tive property allows the RSNS algorithm to be implemented into existing antenna arrays 


and only requires modifying the antenna signal processor. 


Computer simulations of a two-channel RSNS-based DF interferometer have 
shown the effectiveness of the variable resolution signal processing algorithm. That is, 
the multiple levels and the progressive angle strategy can overcome the low SNR restric- 


tions with single element spacing. 
Suggested future work includes: 


e Simulating and conducting the variable analog-to-digital function in order to 


employ the variable resolution signal processing. 


e Implementing the variable resolution signal processing algorithm into a firm- 
ware architecture to achieve the objective of a digital RSNS-based DF sys- 


tem. 


e Developing a more sophisticated noise model, such as mutual coupling inter- 
ference, into RSNS fundamental equations in order to give a more realistic 
signal processing simulation. 
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APPENDIX: SIMULATION CODES 


The programs in this thesis are written in the MATLAB software. All program 


files are listed and briefly explained as follows: 


li. 


findDynamicRange.m is used to find the dynamic range and the discrete states vec- 


tors. 

findPhi.m is used to find the phase adjustment of the modular channel. 
enCode.m is used to quantize the modular folding waveforms into discrete states. 
moduli.m is used to compute the estimated AOA. 

estCorBin.m is used to compute the probability of correct bin estimation. 


VResolution.m is used to demonstrate the three-level variable resolution signal 


processing algorithm. 
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ole 





oe 


= Program Name: findDynamicRnage.m 





oe 





ole 


function [RSNS]= findDynamicRange (modulus) 
[mR mC] =size (modulus) ; 
m=sort (modulus) ; 





cN=mC; % cN=channelNumber 
for i=1:cN 
s(i)=i-1; % s=shift 
end 
period=1; 
for i=1:cN 
period=LCM (period,m(i)); % period=LCM (2*ChannelNumber*modulus) 
end 


period=2*cN*period; 


for i=1:cN 
for j=l+s(i) : cN*m(i)+s(i)+2 


OM(i,j)=floor( (j -s(i)-1)/cN ); % dM=discrete Modulus 


end 

for j=cN*m(i)+3+s(i) : 2*cN*m(i)+s(i) 
AM(i,j)=floor( (2*cN*m(i)+cN-j+s(i))/cN ); 

end 

k=1; 

index=1; 

while (index < 1.5*period) 


OM( i,2*k*cN*m(i)+S(i)+1 : 2* (k+1)*cN¥*m(i)+s(i) )=... 


AM( i,l+s(i) : 2*cN*m(i)+s(i) ); 
index=2* (k+1) *cN*m(i)+s (i); 
k=k+4+1; 
end 
end 
aM=dM' ; 
gb(:, [2:cN+1])=dM(:, [1:cN]) ; 
[sgbr sgbc]=size(gb) ; 
gb(:,1)=(0:1:sgbr-1)'; 
gc=gb; 
t = clock; % calculate processing time. 
max=0; 
maxMatrix=[]; 
k=1; 
while 1 % negelect all zero component 
if gc(k, [2:cN+1] )==zeros(1,cN) 
k=k+4+1; 
else 
break 


j=k+1; 


for i=k:period % calculate dynamic range. ---- 


while j< 1.5*period 
dylId=j; 
flag=0; 
while dyId >i 
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58. if gc(j, [2:cN+1]) == gc(dyId-1, [2:cN+1] ) 


59. flag=1; 
60. break 
61. end 

62. dyId=dyId-1; 
63. end 

64. if flag 

65. break 

66. else 

67. j=j+1; 

68. end 

69. end 

70. h(i,1)=gce(i+1,1); 
71. h(i,2)=9c(j,1); 
72. end 


73. max=0; 
74. maxMatrix=[]; 
75. for i=k:period 


76. dy=h(i,2)-h(i,1)+1; 

77. if dy > max 

78. max=dy ; 

79. maxMatrix=gc(h(i,1):h(i,2), [1:cN+1]); 
80. end 

81. end 


82. timeRead=etime(clock,t) ; 

83. gd(:, [2:cN+2] )=maxMatrix(:, [1:cN+1]) ; 
84. gd(:,1)=(0:1:max-1)'; 

85. gd(max+1,1)=max; 

86. gd(max+1,2)=period; 

87. for k=3:cN+2 


88. gd (max+1,k)=m(k-2) ; 

89. end 

90. RSNS=gd; 

OL. eos eee SSS ete See Ste End of Program------------------------------- 
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= Program Name: findPhi.m 





oe 





oe 


function [phi] =findPhi (binMatrix) 
[bMr bMc] =size (binMatrix) ; 


CN =bMc-2; 
M=binMatrix(bMr,1); 
for i=1:cN 


Schannel 


Number 


sdynamic range value 


m(i)=binMatrix(bMr,i+2); 


end 

index=ceil( M/2 ); 

if M <= index+cN 
setPts=M; 
setLth=M-index; 

else 
set Pts=index+cN; 
setLth=cN; 

end 

for i=1:cN 





table=modulusTable(m(i),cN) ; 
temp=binMatrix (index:setPts,i+2); 


temp=temp!'; 


for j=1:2*cN*m(i) 


comp=table(1 


,jij+setLth) ; 


if temp == comp 


if( rem(M, 


2) ) 


phi (i)=table(2,j); 


else 





phi (i) =(table(2,j)+table(2,3+1))/2; 


end 
break 


60 





ole 








ole 


= Program Name: enCode.m 


oe 








ole 


function [bin] = enCode (Modulus,MixerOp) 
for i=1:Modulus 

ec m(i)=0; 
end 


for i=1:Modulus 


threshold (i) =cos ((2*Modulus-2*1i+1) *pi/ (2*Modulus) ) ; 


end 
for j=1:Modulus 








if ( MixerOp > threshold(j) 


c_m(j) = 1; 
else 


end 
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52. 
53. 





ole 
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= Program Name: moduli.m 





oe 





ole 


function [eAOA] =moduli (ModuliiInf,angle, Limit) 
c=3e8; 

£=8000e6; sradio frequency. 
lambda=c/f; $lambda == wavelength 
k=2*pi/lambda; 
rad=pi/180; 
errorCode=-1; 
epsilo=0.03; 
alpha=0.9; 
mm=Modulilnf (1 
sf=Modulilnf (3 
snr=Modulilnf ( 
theta=angle; 
thetaR=theta*rad; 


S23 


); 
4); 





binMatrix=findDynamicRange (mm) ; 
[bMr bMc] =size (binMatrix) ; 


M=binMatrix(bMr,1) ; sdynamic range. 

cN=bMc-2; channel number. 

nf=(M/ (2*cN) ) ./mm; nf == number of folds 

ref18) 

d=nf.* (lambda/2)*sf; $d == element spacing distance (scaled) 


ref18) 
xi=findPhi(binMatrix); %Phase Adjustment. 
bin=(binMatrix(1:M,1))'; 
%* Assume the reference antenna phase p_ r=0 

% that is, V_I(p_r)=1 & V_Q(p _r)=0 

% then Eq(4.14) V_mI (theta) =cos (k*d_i*sin(theta) )=V_I 
% Eq(4.15) V_mQ(theta)=sin(k*d_i*sin(theta) ) = ( 


if (sf <1) 

Fov=[-90 90]; 
else 
Fov=asin(1/sf) /rad; 
Fov=[(-Fov+1) (Fov-1)]; 





end 

if (Limit(2) <= Fov(1) || Limit(1) >= Fov(2) ) 
eAOA (1:4) =0; 
return; 

else 


Limit (1) =max (Limit (1),Fov(1))j; 

Limit (2) =min(Limit (2) ,Fov(2)); 

Limit=Limit*rad; 
end 
for i=1:cN 

V_Limit (1,1) =cos ( k*d(i)*sin(Limit (1) ) )*cos (x1i(i)) 
k*d(i)*sin(Limit(1)) )*sin(xi(i)); 

V_Limit (1,2) =cos ( k*d(i)*sin (Limit (2) ) ) *cos (x1i(i)) 
k*d(i)*sin(Limit(2)) )*sin(xi(i)); 

env (1,i) =enCode (mm(i),V_ Limit (i,1)); 

env (2,i) =enCode (mm(i),V_ Limit (i,2)); 
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sscale factor. (change here) 


(12- 


(13> 


54. end 
55. for i=1:2 





56. flag=1; 

ow ae for j=1:M 

58. if( env(i,:) == binMatrix(j,3:2+cN) ) 

59. binLimit (i) =binMatrix(j,1) ; 

60. flag=0; 

61. break; 

62. end 

63. end 

64. if flag 

65. bin=errorCode; 

66. disp('out of FOV'); 

67. return; 

68. end 

69. end 

70. %$--------------- Noise Free-------------------------- 

71. for i=1:cN 

72. V_I (i) =cos (k*d(i)*sin(thetaR) ) ; % I term of channel Noise Free 

73. V_Q(i) =sin(k*d(i)*sin(thetaR) ) ; % Q term of channel Noise Free 

74. end 

75. $--------------------------- 5-5-5 ee ee ---- -- ---------e 

76. if (sf < 1) 

77. sf=1; 

78. end 

79. N=100; $number of trials 

80. %--------------- Additive Noise ---------------------- 

81. while (1) 

82. for i=1:cN 

83. VI(i,:)=V_I(i) *ones (1,N) ; 

84. VO(i,:)=V_Q(i) *ones (1,N) ; 

85. end 

86. comVl=complex(VI(1,:),VQ(1,:)); 

87. comV2=complex(VI(2,:),VQ(2,:)); 

88. NcomV1 = awgn(comV1,snr) ; % noise added complex phasor 
for channel 1 

89. Ncomv2 = awgn(comV2,snr) ; 

90. NV(1,:) =cos (xi (1) )*real (NcomV1) -sin(xi(1))*imag(NcomvV1) ; 

91. NV(2,:) =cos (xi (2) ) *real (NcomV2) -sin(xi(2)) *imag(NcomvV2) ; 

92. for n=1:N 

93. NenV(1,n) =enCode (mm(1),NV(1,n)); 

94. NenvV (2,n) =enCode (mm(2) ,NV(2,n)); 

95. end 

96. enCmx=Nenv' ; 

97. for i=1:N 

98. flag=1; 

99. for j=1:M 

100. if ( enCmx(i,1:cN) == binMatrix(j,3:2+cN) ) 

101. Nbin (i) =binMatrix(j,1); 

102. flag=0; 

103. break; 

104. end 

105. end 

106. if flag 

107. Nbin(i) =errorCode; 

108. end 

109. end 
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110. $----------------------------------- 5 --- -- 5-5 5-5 

111. %$----------- probalility ----------------------------- 

112. for i=1:M 

T3. hit (i) =length( find(Nbin==bin(i)) ); 

114. end 

115. hit (M+1) -=length( find(Nbin==errorCode) ); 

116. [num posBin] =max (hit) ; 

117: prob=num/N; 

118. newN=floor( prob* (1-prob) / ((1-alpha) * (epsilo*2))); 

119. if newN >N 

120. N=newN; 

121. Vi=e[]; 

122. VO=[]; 

123. comvVi1=[]; 

124. comv2=[]; 

125. NcomvVi1=[]; 

126. Ncomv2=[]; 

127. NV=[]; 

128. Nenv=[] ; 

129. enCmx= [] ; 

130. Nbin=[]; 

13:1. 3 else 

132. Mb=posBin-1; 

133. if ( binLimit(1) <= Mb && Mb <= binbLimit(2) ) 

134. eAOA (1) =1; % ID Bin 

135. eAOA(2)=asin( (2*Mb+1) / (M*sf)-1/sf£) /rad; % Main AOA 

136. binWidthe= ( asin ((2*Mb-M+2) / (M*sf) ) - asin ( (2*Mb- 
M)/(M*sf)) )/rad; 

137". eAOA (3) = eAOA(2) -binWidth/2; % NextLim(1) 

138. eAOA (4) = eAOA(2)+binWidth/2; % NextLim(2) 

139. else 

140. eAOA (1:4) =0; 

141. end 

142. break; 

143. end 

144. end 

145. GeSSsssSese Seer eeeS End. of -Program=-=-=s5+s-SsssSs5sss3SSse5Ss555 
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= Program Name: estCorBin.m 





oe 





ole 


function [Prob] =estCorBin (ModulilInf, angle) 


c=3e8; 

£=8000e6; sradio frequency. 
lambda=c/f; Slambda == wavelength 
k=2*pi/lambda; 

rad=pi/180; 


errorCode=-1; 
epsilo=0.03; 
alpha=0.9; 
mm=Modulilinf (1 
sf=Modulilnf (3 
snr=Modulilnf ( 
theta=angle; 
thetaR=theta*rad; 


SZ 
y $scale factor. (change here) 
4); 





binMatrix=findDynamicRange (mm) ; 
[bMr bMc] =size (binMatrix) ; 


M=binMatrix(bMr,1) ; sdynamic range. 
cN=bMc-2; channel number. 
nf=(M/ (2*cN) ) ./mm; snf == number of folds (12- 
ref18) 
d=nf.* (lambda/2)*sf; $d == element spacing distance(scaled) (13- 
ref18) 
xXi=findPhi(binMatrix); %Phase Adjustment. 
bin=(binMatrix(1:M,1))'; 
%* Assume the reference antenna phase p_r=0 
% that is, V_I(p_r)=1 & V_Q(p_r)=0 
% then Eq(4.14) V_mI (theta) =cos(k*d_i*sin(theta) )=V_I(p_m) 
% Eq(4.15) V_mQ(theta) =sin(k*d_i*sin(theta) )=V_Q(p_m) 
Sars eH SSeS SSS Noise Free-------------------------- 
Tbin=[]; 
for i=1:cN 
V_I (i) =cos (k*d(i)*sin(thetaR) ) ; % I term of channel Noise Free 
V_Q(i) =sin(k*d(i)*sin(thetaR) ) ; % Q term of channel Noise Free 
V(i) =V_I (i) *cos (xi (i))-V_Q(1i) *sin(xi(i)); 
env (1) =enCode (mm(i),V(i)) ; 
end 
for j=1:M 
if( env(1:2) == binMatrix(j,3:2+cN) ) 
Tbhin=binMatrix(j,1); 
break; 
end 
end 
if (sf «< 1) 
sf=1; 
end 
N=100; $number of trials 
$--------------- Additive Noise ---------------------- 
while (1) 


for i=1:cN 
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54. 
55:2 
56. 
oa 
58. 
59. 


60. 
61. 
62. 
63. 
64. 
65. 
66. 
67. 
68. 
69. 
70. 
TA 
72. 
73. 
74. 
75. 
76. 
UT 
78. 
79. 
80. 
81. 
82. 
83. 
84. 
85's 
86. 
87. 
88. 
89. 
90. 
91. 
92. 
93. 
94. 
95. 
96. 
97. 
98. 
99. 


100. 
OA. 
102. 
103. 


VI(i,:)=V_I(i) *ones (1,N) ; 
VQ(i,:)=V_Q(i) *ones(1,N) ; 
end 
comVl=complex(VI(1,:),VQ(1,:)); 
comV2=complex(VI(2,:),VQ(2,:)); 
NcomV1 = awgn(comV1,snr) ; % noise added complex phasor 
for channel 1 
NcomV2 = awgn(comv2,snr) ; 
NV(1,:) =cos (xi(1))*real (NcomV1) -sin(xi(1))*imag(NcomvV1) ; 
NV(2,:) =cos (xi (2) ) *real (NcomV2) -sin(xi(2)) *imag(Ncomv2) ; 
for n=1:N 
Nenv(1,n) =enCode (mm(1),NV(1,n)); 
Nenv (2,n) =enCode (mm(2) ,NV(2,n)); 
end 
enCmx=Nenv' ; 
for i=1:N 
flag=1; 
for j=1:M 
if ( enCmx(i,1:cN) == binMatrix(j,3:2+cN) ) 
Nbin (i) =binMatrix(j,1); 
flag=0; 
break; 
end 
end 
if flag 
Nbin (i) =errorCode; 
end 


a probalility ----------------------------- 
hit=length( find(Nbin==Tbin) ); 
prob=hit/N; 
newN=floor( prob* (1-prob) / ((1-alpha) * (epsilo*2))); 
if newN > N 

N=newN; 

VI=[]; 

VO=[]; 

comV1=[]; 

comv2=[]; 

NcomvV1=[]; 

Ncomv2=[]; 

NV=[] ; 

Nenv=[]; 

enCmx= [] ; 

Nbin=[]; 
else 

Prob=prob*100; 

break; 
end 
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= Program Name: VResolution.m 





oe 





oe 


cle 
clear 
fid = fopen('Figure37v0.txt','w'); 
modu0309= [3 9 105/112 20]; 
modu0816=[8 16 1 20]; 
modu0824=[8 24 105/102 20]; 
thetaMatrix=-70:1:70; 
for i=1:length(thetaMatrix) 
angle=thetaMatrix (i) 
Limit=[-72 72]; 
Vi (243 ):=0% 
eAOA0309=moduli (modu0309,angle, Limit) ; 
eAOAref =moduli(modu0816,angle,Limit) ; 
y (1) =angle-eAOA0309 (2) ; 
yRef=angle- eAOAref (2) ; 
if (eAOA0309 (1) ==1) 
Limit=eA0A0309 (3:4); 
e€AOA0816=moduli (modu0816,angle, Limit) ; 
y (2) =angle-eAOA0816 (2) ; 
if (eAOA0816 (1) ==0) 


y(2)=y(1); 
y(3)=y(1); 
else 


Limit=eAOA0816 (3:4) ; 
eAOA0824=moduli (modu0824,angle, Limit) ; 
y (3) =angle-eAOA0824 (2) ; 
if (eAOA08 24 (1) ==0) 
y(3)=y(2); 
end 
end 
end 


fprintf (fid, '$d\t%5.2£\t%5.2£\t$5.2£\ts5.2£\n',angle,y(1:3) ,yRef) ; 
end 
fclose(fid); 


g 


fid = fopen('Figure37v0.txt'); 
a = fscanf (fid, 'sd\tsf\ts£F\cts£F\tst', [5 inf]); 





a vas as 

fclose(fid) ; 

[Nr Nc] =size(a); 

x =( a( 1:Nr,1) )'; 

yl =( a( 1:Nr,2) )'; 

y2 =( a( 1:Nr,3) )'; 

y3 =( a( 1:Nr,4) )'; 

y4 =( a( 1:Nr,5) )'; 

figure (1) 

plot (x,y4, 'LineWidth',2); 

grid; 

xlabel('Angle of Arrival (degrees)','FontSize',14); 
ylabel('Estimative Error (degrees)','FontSize',14); 


set (gca,'YTick',-60:10:60, 'FontSize',14) ; 
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print -f1 -tiff -depsc figure37v0 
figure (2) 

plot (x,y1,'LineWidth',2); 
xlabel('Angle of Arrival (degrees) ' 
ylabel('Estimative Error (degrees) ' 
set (gca, 'YTick',-15:5:15,'FontSize' 
grid; 

print -f2 -tiff -depsc figure38v0 
figure (3) 

subplot (2,1,1); 

plot (x,y1,'LineWidth',2); 
xlabel('Angle of Arrival (degrees) ' 
ylabel('Estimative Error (degrees) ' 
set (gca, 'YTick',-15:5:15,'FontSize' 
grid; 

subplot (2,1,2); 

plot (x,y2, 'LineWidth',2); 

grid; 

xlabel('Angle of Arrival (degrees) ' 
ylabel('Estimative Error (degrees) ' 
set (gca, 'YTick',-15:5:15,'FontSize' 
print -f3 -tiff -depsc figure39v0 
figure (4) 

subplot (2,1,1); 

plot (x,y2, 'LineWidth',2); 
xlabel('Angle of Arrival (degrees) ' 
ylabel('Estimative Error (degrees) ' 
set (gca, 'YTick',-15:5:15, 'FontSize' 
grid; 

subplot (2,1,2); 

plot (x,y3,'LineWidth', 2); 

grid; 

xlabel('Angle of Arrival (degrees) ' 
ylabel('Estimative Error (degrees) ' 
set (gca, 'YTick',-15:5:15,'FontSize' 
print -f4 -tiff -depsc figure40v0 
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,'FontSize'! 
,'FontSize'! 
,14); 


,'FontSize'! 
,'FontSize'! 
,14); 


,'FontSize'! 
,'FontSize'! 
, 14); 


,'FontSize'! 
,'FontSize' 
,14); 


,'FontSize'! 
,'FontSize'! 
, 14); 


ISS SSeS Sse SSeS End of Program-------------------------- rrr 
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